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ABSTRACT 

We investigate the behaviour of asymptotic giant branch (AGB) stars between metal- 
i-S^ ■ hcities Z = 10"'' and Z = 10~®. We determine which stars undergo an episode of 

flash-driven mixing, where protons are ingested into the intershell convection zone, 
Q , as they enter the thermally pulsing AGB phase and which undergo third dredge-up. 

^—^ • We find that flash-driven mixing does not occur above a metallicity of Z = 10^^ for 

jj^ I any mass of star and that stars above 2 Mq do not experience this phenomenon at 

d . any metallicity. We find carbon ingestion (CI), the mixing of carbon into the tail of 

hydrogen burning region, occurs in the mass range 2 M© to around 4Mq. We suggest 
that CI may be a weak version of the flash-driven mechanism. We also investigate 
^ . the effects of convective overshooting on the behaviour of these objects. Our models 

' struggle to explain the frequency of CEMP stars that have both significant carbon 

fSJ , and nitrogen enhancement. Carbon can be enhanced through flash-driven mixing, CI 

CO ' or just third dredge up. Nitrogen can be enhanced through hot bottom burning and 

' the occurrence of hot dredge-up also converts carbon into nitrogen. The C/N ratio 

. may be a good indicator of the mass of the primary AGB stars. 
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^ . 1 INTRODUCTION nary systems l|Lucatello et al.ll2005l ). This is consistent with 

5-H ' , , , . , , , . . , the whole population being binaries. 
' In recent years, there has been considerable interest m the 

■ ■ ■ ' observation of metal-poor stars and particularly those of a This scenario may be compUcated somewhat. It has 
carbon-rich nature. A surprising fraction of stars with metal- commonly been assumed that, when the material is accreted 
licities below [Fe/H] of around -2.5 are found to be rich in on to the companion, it remains on t h e surf ace of the re- 
carbon. This fractio n, while uncertain, m ay be as high as cipient star. However, IStancliffe et^ (|2007|) showed that 
around 20 per cent (|Lucatello et al.ll2006l ). Many of these thermohaline mixing can lead to this accreted material be- 
stars display evidence of s-process enrichments. It has been ing extensively mixed with the pristine matter of the com- 
proposed that they may be formed by mass transfer in bi- panion and this can have important consequences for the 
nary systems. Such a system would consists of two stars in abundances that would be observed in such an object. The 
a wide orbit. The primary star evolves to the asymptotic extent of this mixing depends on the composition and mass 
giant branch (AGB) where it becomes enriched in carbon of accreted mater ial, as well as the evolutiona ry state of the 
(and also s-process elements). Strong mass loss strips the accreting star (seelStandiffe fc Glebbee3l2008|). However, we 
envelope from this star which becomes a white dwarf and may still hope to observe the signatures of low-metallicity 
fades from view. The star is losing mass so some of this is AGB stars in the compositions of the CEMP stars that are 
accreted by the companion which becomes carbon-rich. It is still visible today. To do this demands that we know what 
this companion that we see today. Credence is lent to this ar- the expected nucleosynthetic signatures of the parent AGB 
gument by the observation that, of the carbon and s-process stars are. 

rich stars, around 70 per cent are found to be members of bi- The asymptotic giant branch (AGB) is a late phase 

of stellar evolution for stars with masses from 1 to around 
SMq. It follows the end of the core helium burning phase. 

* hlau@kiaa.pku.edu.cn As the core of the star runs out of helium it begins to ascend 
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the AGB. The star first expands and cools and the convec- 
tive envelope deepens. The envelope can reach down as far 
as the hydrogen burning shell, pulling material to the sur- 
face and altering the surface composition of the star. This 
is the second dredge-up 

Helium and hydrogen are both burning in thin shells. 
The helium burning shell slowly moves outwards in mass 
and the region between the two shells narrows. When they 
get very close (separated by a few hundredths of a solar 
meiss) the star enters the thermally pulsing asymptotic gi- 
ant branch (TP-AGB) phase. Thermal pulses are the rapid 
increase in the helium-burning luminosity to over 10* Lq , 
for a brief period of time, usually around lOyr. The ther- 
mal pulses are a consequence of the thinness of the he- 
lium burning shell and high temperature sensitivity of the 
tripl e-g reaction. The theory is described, for AGB stars, by 
ISchw arzschild & Harm (1965i). When energy is dumped in 
the thin shell, the temperature can rise when the shell ex- 
pands. There is no thermostatic control until the shell is no 
longer thin. Also, for the pulse to occur, the increased radia- 
tive losses owing to the increase of shell temperature must 
not carry energy away faster than it is being generated. This 
is usually fulfille d because helium burning is highly sen sitive 
to temperature. lYoon. Langer fc van der Sluvd l|2004) have 
investigated the stability of shell sources and find that the 
helium burning shell tends to be more stable at high tem- 
perature because the triple-a reaction is less temperature 
sensitive. 

During the thermal pulses, the sudden increase in en- 
ergy release, caused by enhanced helium burning, leads to 
convection between the two burning shells and mixes the 
ashes of helium burning outwards from the helium burning 
shell. The star expands and the hydrogen burning shell cools 
and may be extinguished. Eventually thermostatic control is 
restored in the He-shell, the helium luminosity starts to drop 
and the intershell convection region shuts down. The convec- 
tive envelope deepens during this power down phase, which 
lasts for about 100 yr. The envelope could penetrate into re- 

5 ions wher e the intershell convection zone had been active 
IbenlllQTSi ). This means that the carbon that has been pro- 
duced during helium burning is ingested by the envelope. 
This is called third dredge-up (TDUP). The surface is then 
enriched in carbon (and s-process element^. The base of 
the convective envelope can be hot enough for some CNO 
cycle reactions to occur, particularly for stars more massive 
than AMq. This is called hot bottom burning (HBB). He- 
lium burning returns to a quiescent state and the hydrogen 
shell reignites. Helium accumulates again during this inter- 
pulse period which can last over 1,000 yr. The cycle then 
repeats and may do so many times before stellar winds can 
strip away the envelope. 

Low-metallicity AGB stars are known to evolve dif- 
ferently to their higher metallicity counterparts. In the 
case of the low-mass objects it is possible for the helium- 
burning driven intershell c onvection zone to pen etrate into 
the hydrogen-rich envelope (|Fuiimoto et al.|[2000l ). This does 



not happen at higher metallicity because the hydrogen burn- 
ing shell presents an effective entropy barrier which inhibits 
the advance of the intershell convection zone. As the metal 
content of the star is reduced, hydrogen burning via the 
CNO cycle is weaker and presents less of a barrier. Be- 
low some critical metallicity, the intershell zone manages to 
penetrate into H-rich regions. This leads to hydrogen being 
dragged down to regions of high temperature, with vigorous 
H-burning occurring as a result. The convective region can 
then split in to two separate zones, one driven by H-burning, 
the other by He-burning. This phenomenon is referred to as 
fiash driven mixing (FDM, ISuda et al.ll2004l). It has been le- 
por ted by several authors, i nclud ing Fuiimoto et al.l (|2000l ) 
and lCampbeU fc Lattanriol (|2008l )P[ 

Metal-free evolution for lo w- and inter m ediate -mass 
stars has been explored b y 'Chief fi et al.l (1200 ll ) and 
Sies s et al.l (|2002l '). ICliieffi et a l. (2001) find that during ther- 
mal pulses of AGB stars at or below 5 M0 , a convective shell 
forms at the H-He interface and the subsequent expansion 
of this convective shell in to the underlying layers previously 
occupied by the pulse dredges up some carbon and initiates 
an H fiash. In the p aper, this phenom enon is referred as car- 
bon mg estion (Cl). ISiess etaU l|2002h also find this unusual 
convective shell and that there is a mixing episode for stars 
at or below 6 Mp . They a lso re fer to it as CI. 

In addition, iHerwid l|2004 ) reported that, at a metal- 
licity oi Z = 10~*, the nature of third dredge-up changes. 
With the inclusion of convective overshooting, he found that 
hydrogen would actually burn during TDUP and that the 
propagation of a H-burning fiame could erode the entire in- 
tershell and even dredge up some material from the CO core. 
This behaviour has been dubbed hot dredge-up. 

Our purpose in this work is to delineate the boundaries 
between these regimes of behaviour. This will enable us to 
place limits on which stars are able to produce which nucle- 
osynthetic signatures. Based on our models, we can hope to 
explain the variation of carbon and nitrogen enrichment of 
CEMP stars in terms of the different enrichment processes of 
extremely metal-poor AGB stars of different masses. There- 
fore, the observed frequencies of different types of CEMP 
stars may tell us the relative numbers of low- and high-mass 
AGB stars in extremely metal-poor environments. This can 
help deduce the fo rm of the initial ma ss function at low 
metallicity (see e.g. iLucatello et al.ll2005h . 



2 STELLAR EVOLUTION CODE 

We use the stellar evolution code S TARS whic h is a vari- 
ant of that originally develope d by Eggleton (|l97ll ) and 
updated by many authors (e.g. iPols et al.lll995h . The ver- 
sion used here includes the A GB-specific modifications of 
IStancliffe, Tout fc Polsl (|2004|). together with the updated 
opacitv tables of lEldridge fc Tout! l|2004 ) which account for 
the enhancement to the opacity from enhancements in car- 
bon and oxygen abundance. The code has been extensively 



^ The s-process is the capture of neutrons by nuclei in conditions 
of sufficiently low neutron density that any unstable products 
have time to decay before another neutron is captured. The neu- 
tron source in AGB stars is not yet properly understood. 



^ The community is yet to accept a standard nomenclature 
for these event s. They have also b een called 'proton inges- 
tion episodes' jCrista llo et al.l |2007|) and 'dual shell flashes' 
llCampbell fc Lattanzidl2008l) . 
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Figure 1. AGB behaviour of tlie models as a function of mass and metallicity witliout oversliooting. The symbols are T — the model 
undergoes third dredge-up, H - the model undergoes hot third dredge-up, F — the model undergoes flash-driven mixing, CI - the model 
undergoes a carbon injection, X — no TDUP occurs, S - the model undergoes carbon ignition before the AGB and presumably becomes 
an S-AGB star. 



used to compute TP- AGB evolution i n recent years (see e.g. 
IStancliffe et alll2005l : iLau et al.|[2007l l 

We do not apply mass loss at any stage of evolution. 
Mass-loss rates for AGB stars are very uncertain and it is ex- 
pected that mass loss scales with metallicity. In any case, the 
behaviour we are interested in happens in the first few pulses 
before any mass loss would be significan t. We use a = 2.0 
for the mixing length parameter in the Bohm-Vitens3 (|l958l ) 
prescription. This value is chosen based on calibration to a 
solar model. We also cho se solar-scaled initial abundances 
ijAnders fc Gre vessel [l989l ) rather than attempting to follow 
the chemical evolution of the Galaxy. 

We have made models of metallicity Z = 10~'*,3 x 
10"^10~^10"^10"■',10"** and 0. We have not modelled 
the lowest-mass stars at the lowest metallicities because 
these stars undergo a helium flash. The STARS code is cur- 
rently unsuitable for modelling this phase in a star's evolu- 
tion . While it is possible to g enerate a post-He flash model 
(see IStancliffe &i Jeffervll2007l . for details), the method as- 
sumes that the envelope composition is unaffected by the 



helium flash and this c annot be guaranteed below a metal- 
licity of [Fe/H] < -4.5 ijFuiimoto et al.lll990l l. 

The models are evolved from the pre-main sequence to 
the TP-AGB with 999 mesh points. At the beginning of 
the TP-AGB, the AG B-specific mesh spacing function of 
IStancliffe etahl l|2004l ) is switched on. It is also at this point 
that the AGB mixing prescription by the same authors is 
applied. In the initial runs, no convective overshooting was 
applied at any phase in the star's evolution. 



3 RESULTS 

We can roughly divide the evolutionary behaviour of these 
extremely metal-poor stars into five different regimes, (i) 
stars with a fiash-driven mixing episode, (ii) stars with CI, 
(iii) stars with no third dredge-up, (iv) stars that experi- 
ence third dredge-up which can be subdivided depending on 
whether hot dredge-up or hot bottom burning or both are 
present and (v) stars which undergo carbon ignition before 
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Figure 2. The abundance of CNO elements by mass at the beginning of the TP-AGB. The action of second dredge-up substantially 
raises the abundance in the case of the intermediate-mass stars. 



0.61 



0.6 



0.59 



0.58 



0.57 



H-burning shell 



5x10" 



Convective region 



He-burning shell 



C/0 core 



5.5x10" 6x10" 
Model Number 




6.5x10" 



Figure 3. Evolution of 1.5 Mq, Z = 10^^ model during He-FDM. We can see the convective region reaches the H-burning shell at model 
number around 56000. The messy appearance is due to numerical problems to resolve the mixing episode. The 2nd pulse, at about model 
62,000 is simply a normal thermal pulse. 
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the AGB. These regimes of behaviour are shown in Fig.[T]in 
msiss and metallicity space. 

Stars with an initial mass of 7Mq or above be- 
come super-AGB (SAGB) when carbon i s ignited non- 
degen erately before any thermal pulses l|Gil-Pons et al.l 
I2OO5I ). We do not follow the evolution of such SAGB stars, 
owing to t he consider able complexity involved in their com- 
putation (|Siessll2006l ). It is worth noting that a different 
mixing-length parameter a could change the critical mass at 
which a star becomes an SAGB star. For example, a 7Mq 
model oi Z = 10~* evolved with a = 1.925 does not evolve 
into a SAGB star. This is a perfect example to show that 
the mixing length parameter is crucial to determine the evo- 
lution of extremely-metal poor stars (and in fact, any star). 
In this work all models are computed with a — 2.0 and we 
shall examine the effect of varying it in the future. 

Stars with Z = 10~* (corresponding to [Fe/R] ^ -2.3) 
behave similarly to their higher metallicity cousins, with the 
occurrence of third dredge up and hot bottom burning. This 
suggests that AGB stars with Z — 10"* have enough CNO 
elements in the hydrogen burning shell that their evolution 
is not affected by the low metallicity. When the metallic- 
ity is decreased to Z — 10^^ — 10"'', the high mass AGB 
stars still behave as they do at higher metallicity, in that 
they experience only thermal pulses, third dredge-up and 
hot bottom burning. This is because second dredge-up is 
more efficient in higher mass stars and leads to an increased 
CNO abundance in the envelope and (more importantly) the 
hydrogen burning shell (see Fig. [2}. Therefore, the amount 
of metals, particularly the CNO elements in the hydrogen 
burning shell is one of the important factors in determining 
the evolutionary behaviour. 



3.1 Helium-flash-driven mixing and carbon 
ingestion 

For stars with metallicity at or belo w 10~^ and masses a t 
or below 2M0 He-FDM can occur l|Fuiimoto et al.lll990l l. 
He-FDM occurs when a convective zone, driven by helium 
burning, penetrates into hydrogen-rich regions mixing pro- 
tons into the carbon-rich parts of the star. This can happen 
either during the helium flash on the red giant branch or 
during the TP-AGB phase. It is the latter case that con- 
cerns us here. The outer edge of the intershell convective 
zone (ICZ), generated by helium burning, extends into lay- 
ers containing hydrogen. Hydrogen is then ingested by the 
convective zone and is mixed down to a depth where the tem- 
perature is much hotter than the hydrogen burning shell. As 
a result, the hydrogen-burning luminosity increases substan- 
tially. This leads to the convective zone splitting into two 
parts. The upper convective shell is now driven by hydrogen 
burning, while the lower convective shell is still driven by 
helium burning. As more hydrogen is mixed into the hydro- 
gen convective shell from above, the hydrogen flash develops 
further. After the decay of the flash, the region formerly oc- 
cupied by the hydrogen-driven convective shell is carbon- 
and nitrogen-enhanced because of the mixing of material 
from the helium-burning shell. Then the base of the convec- 
tive envelope penetrates into the layers formerly occupied by 
the hydrogen burning shell and brings carbon and nitrogen 
to the surface. 

Fig. [3] shows an example of a star undergoing FDM, in 
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Figure 4. Top panel: helium luminosity as a function of time 
since the first thermal pulse. Middle panel: H-exhausted core mass 
(dotted line) and He-exhausted core mass (solid line) line as a 
function of time. Bottom panel: abundance by mass fraction of 
carbon (solid line), nitrogen (dotted line) and oxygen (dashed 
line). These plots are for the 1.5 model at Z = 10"^. The 
model undergoes an episode of FDM after about 2.5 X 10^ yr 
since the first thermal pulse, leading to enhancements in the CNO 
abundances. Thereafter, only carbon and oxygen are enhanced. 



this case a 1.5 Mq star of metallicity Z = 10~^. At around 
model number 55000, the convective zone, driven by helium 
burning, penetrates the hydrogen-burning shell. This causes 
the hydrogen luminosity to peak. The increased hydrogen 
luminosity drives a convective zone which connects with the 
intershell convective zone. Material is then mixed between 
the hydrogen burning shell and the top of helium burn- 
ing shells and so there is significant carbon present in the 
hydrogen-burning shell. Later, around model number 60000, 
the convective envelope deepens and reaches the region that 
has been previously mixed. In this way, the products of he- 
lium burning are dredged up to the surface during the He- 
FDM episode. Nitrogen is also dredged up because the CNO 
cycle converts carbon into nitrogen. 

Fig. |4] shows that the surface abundances of CNO 
elements are enhanced significantly during the He-FDM 
episode, corresponding to the third pulse of the 1.5 Mq, 
Z — 10~^ model. The surface abundance of carbon increases 
by more than 3 orders of magnitude to the order of 10~^ be- 
cause helium burnt materials are dredged up to the surface. 
Nitrogen increases from about 2 x lO"'^ to 10~* while oxy- 
gen increases by one order of magnitude to about 4 x 10~^. 
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Figure 5. Top panel: helium luminosity as a function of time 
since the first thermal pulse. Middle panel: H-exhausted core mass 
(dotted line) and He-exhausted core mass (solid line) line as a 
function of time. Bottom panel: abundance by mass fraction of 
carbon (solid line), nitrogen (dotted line) and oxygen (dashed 
line). These plots are for the SMq model a,t Z = 10"'^. This 
model suffers the CI mechanism and the first episode of TDUP 
leads to enhancements of C, N and O. 



It is also worth noting that after the He-FDM episode, the 
surface CNO abundances resemble those of a much higher 
metallicity star. As a result of higher CNO abundances in 
the hydrogen burning shell, the helium convective zone can- 
not any longer reach the hydrogen burning shell in the sub- 
sequent pulses and the star behaves similarly to a higher 
metalhcity star. Carbon and oxygen are slowly enhanced 
during subsequent pulses by third dredge up just as in higher 
metallicity stars. 

We find that FDM episodes only occur for the lowest 
stellar masses (2 Mq and below) . Stars of higher mass do not 
have violent enough thermal pulses to drive the ICZ up into 
the H-depleted regions. However they do have pulses that 
are violent enough to trigger TDUP. Once TDUP has raised 
the CNO abundance sufficiently, H-burning in the shell pro- 
ceeds efficiently via the CNO cycle, the H-shell presents an 
effective energy barrier to the ICZ and prevents any FDM. 

At Z = 10~® there is a change in the evolution of the 
3 Mq and 4 Mq star. We find that an extra convective zone, 
located above the hydrogen burning shell, opens after the 
peak helium luminosity of the thermal pulse. Unlike the 
FDM mechanism, this convective zone is not a breakaway 
region from the intershell convection zone, but is distinctly 




Figure 6. FDM hydrogen abundances (solid line) and Vr — Vad 
(dotted line) of our f.5MQ, Z = 10~^ model during flash driven 
mixing. A region is convective if Vr — Vad > 0. The intershell 
convective zone reaches the H-burning shell, where the abundance 
of hydrogen is around O.f . As a result, plenty of protons are mixed 
down to regions with higher temperature and carbon abundances. 
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Figure 7. CI hydrogen abundances (solid line) and Vr — Vad 
(dotted line) of our 3Mq, Z = 1Q~^ model during carbon in- 
gestion. A region is convective if Vr — Vad > 0- The intershell 
convective zone reaches the tail H-burning shell, where the abun- 
dances of hydrogen is around 10~^. Compared to FDM, far fewer 
protons are mixed, though still enough to drive a convective zone 
by the increased hydrogen burning luminosity. 



se parated from i t. Thi s mix ing episode is refe rred to as a CI 
bv lChiefB et"al] ()200l[ ) and lSiess et all ()2002l ). As described 
by these authors, carbon from the intershell is injected into 
the H-burning shell by this convective zone. This mechanism 
also occurs at 3 - 6Mq, Z = 10"'' and 3 - 4Mq, Z = 10"® 
models. An example of a star that undergoes this process is 
shown in Fig. [5] 

We note that there is a distinction between the nature of 
flash-driven mixing (FDM) and carbon injection (CI). The 
FDM phenomenon involves protons being mixed down into 
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Figure 9. Hydrogen luminosity (solid line) and helium luminosity 
(dotted line) evolution of a 1.5 model at Z = 10~^. The 
hydrogen luminosity reaches 10^ Lq during the FDM at around 
model number 56,000. 



the intershell region and burning at very high temperatures. 
For the CI phenomenon, mixing occurs when a convective 
pocket opens up above the H-shell and its lower edge pen- 
etrates down into C-rich regions. H-burning takes place at 
much lower temperatures than FDM. In addition, the FDM 
involves the ingestion of a small quantity of hydrogen into 
a carbon-rich region. Conversely, CI involves the injection 
of a small amount of carbon into an H-rich region while the 
amount of hydrogen mixed is minimal. 

In both cases, a convective zone opens up at the top 
of the hydrogen burning shell, so when the third dredge-up 
deepens, both mechanisms lead to substantial enrichments 
of nitrogen, though both have only a single episode of N- 
enhancement. Once TDUP has occurred the CNO abun- 
dance is high enough that the star evolves like a higher 
metallicity star. Are these two distinct phenomenon? Both 
are triggered by the inges tion of protons in to the intershell 
convection zone (see e.g. ISiess et al.ll2002i . for the case of 
the CI). Perhaps it is just the quantity that is ingested that 
determines which episode a given star experiences. 

We suggest that proton ingestion and He-FDM may 
have the same origin, namely the mixing of protons into a 
carbon rich region but the amount of hydrogen and carbon 
mixed determines whether the two convective zones connect 
each other and hence whether CI or He-FDM occurs. Figs[B] 
and [7] show how the convective zone reaches the hydrogen 
burning shell in the FDM and CI events respectively. The 
convective region penetrates the burning shell and reaches 
regions where the hydrogen abundance is as high as 0.1 dur- 
ing a FDM event, while it only reaches the tail of the burning 
shell, where the hydrogen abundance is about 10~^, in the 
CI episode. In both cases, the mixing of protons and car- 
bon causes an increase of hydrogen luminosity. In the case 
of FDM the burning is much more vigorous so a larger con- 
vective zone is driven and the mixing between the hydrogen 
burning and helium burning layer is more efficient. 

Figs [3] and [8] show the internal structure around hydro- 
gen burning shell of two stars during FDM and CI respec- 
tively. In the case of FDM, the vigorous hydrogen burning 
immediately drives a convective zone that connects with the 
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Figure 10. Hydrogen luminosity (solid line) and helium lumi- 
nosity (dotted line) evolution of a 3Mq model a,t Z = 10~^. The 
hydrogen luminosity rises to a second peak of 10"* Lq (around 
model number 48,000) in the second pulse and this does not oc- 
cur in the other pulses. It corresponds to the CI. 



intershell convective zone. In the case of CI, the burning- 
driven convection zone appears later in time and is never 
connected with intershell convective zone while it is extant. 
FDM mixes material from the H-burning shell with material 
from He-burning shell in a single large convection zone. For 
CI the mixing is a two-step process. As a result, when third 
dredge up occurs later, more carbon is dredged up to the 
surface in a FDM episode than a CI episode. 

In Figs and 1101 we see there is a increase of hydro- 
gen luminosity after a thermal pulse which correspond to 
when FDM and CI occurs. This is unusual. The hydrogen 
luminosity usually decreases monotonically after a thermal 
pulse as in the later pulses. We can see that the increase 
of hydrogen luminosity occurs after the intershell convective 
zone reaches the hydrogen burning shell by comparing the 
model numbers with Figs |3] and [S] The increase of hydrogen 
luminosity starts at model number 55580 for the FDM and 
at model 48672 for the CI. Compared with Figs |3] and H 
we can see that the increases of luminosity just precedes the 
opening of the second convection zone around the hydrogen 
shell. This is because the convective zone around hydrogen 
burning shell is driven by the increased hydrogen luminosity 
following the mixing by intershell convective zone. The in- 
crease of hydrogen luminosity is much higher in FDM than 
CI. The hydrogen luminosity jumps to 10^ Lq in the case 
of FDM and 10^ L© in CI. The peak of the hydrogen lumi- 
nosity is unsurprisingly correlated with the number of pro- 
tons mixed into the burning region. The FDM case shows a 
much stronger peak because the convective shell penetrates 
through the hydrogen burning shell and causes more mixing 
than following the CI episode. 

In terms of the effect on surface abundances, the distinc- 
tion between FDM and CI is not so clear cut in some cases. 
In fact, a model with CI but near the FDM-CI boundary can 
dredge up almost as much carbon as a model with FDM. On 
the other hand, models with CI, but near the boundary with 
normal AGB evolution, dredge no more than in normal third 
dredge-up. In some situations, more nitrogen than carbon is 
dredged up after CI. This is because the star also undergoes 
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Figure 8. Evolution of a 3Mq, Z = 10"'' model during CI. We can see the intershell convective zone reaches the tail H-burning shell 
at model number around 47000. Then a convective region driven by hydrogen burning reaches deep enough to previously mixed region 
by the intershell convective zone around model number 50000. At around model number 50500, the convective envelope deepens and 
connect with the convective zone driven by hydrogen burning. The third dredge up thus dredged up carbon and nitrogen indirectly 



hot bottom burning which convert carbon into nitrogen. In 
general, the effect of CI is stronger when metallicity is lower. 
This could be because at lower metallicity, the abundance 
of carbon at the hydrogen burning shell is lower, so the car- 
bon ingested causes a much larger increase in the hydrogen 
luminosity, which drives a more extensive convective zone. 



3.2 Hot dredge-up 

In stars of 4Mq and above, we find some differences in the 
dredge-up process. For the more massive stars, we note the 
occurrence of what has been referred to as hot third dredge- 
up (HTDUP) . This phenomenon was described in detail by 
iHerwid l|2004l ). Here, the temperature at the base of the con- 
vective envelope at the point when third dredge-up occurs 
is so high that carbon is converted to nitrogen, resulting in 
a simultaneous enhancement of both carbon and nitrogen. 
HTDUP is an extremely efficient dredge-up, with the en- 
velope reaching almost to the bottom of the intershell (or 



perha ps beyond it, if overshooting is included, see iHerwigj 
I2OO4I for further details). The behaviour of an example of 
a star that undergoes this process is shown in Fig. 1111 CI 
occurs at the first episode of third dredge up in this model 
and the CNO abundances are drastically increased. Then 
hot third dredge up sets in after two more pulses, leading to 
simultaneous increase in carbon and nitrogen. Hot bottom 
burning converts carbon into nitrogen and nitrogen becomes 
more abundant than carbon by more than one order of mag- 
nitude. 



3.3 A metallicity threshold for dredge-up? 

As we move to the lowest metallicity, we find that TDUP 
ceases altogether for some models. The 5Mq model, with 
Z = 10~* has very weak thermal pulses and these are not 
strong enough to bring the convective envelope into con- 
tact with the C-rich intershell. On the other hand, the 6 Mq 
model, with Z = 10~^ does have a third dredge up. The 
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Figure 11. Top panel: helium luminosity as a function of time 
since the first thermal pulse. Middle panel: H-exhausted core mass 
(dotted line) and He-exhausted core mass (solid line) line as a 
function of time. Bottom panel: abundance by mass fraction of 
carbon (solid line), nitrogen (dotted line) and oxygen (dashed 
line). These plots are for our 4Mq model at Z = 10~^. The first 
episode of TDUP, while extremely shallow, dramatically enhances 
the CNO abundances. HTDUP sets in after about 25,000 yr, lead- 
ing to a simultaneous increase in carbon and nitrogen before HBB 
takes over in the interpulse. 



CNO content of the 6Mq in the hydrogen burning shell is 
higher than the 5 Mq because of deeper second dredge up. 
This seems to confirm the strength of the pulses are depen- 
dent on the CN O abundan c es in the hydrogen burning shell, 
as suggested bv lLau et al.l l|2008l ). 

The issue of the efficiency of third dredge-up is a peren- 
nial problem for modellers of AGB stars. Different stellar 
evolution codes give very different predictions for the effi- 
ciency of TDUP. Some of the reason for this uncertainty 
is due to the mathematical treatment of the boundaries 
between radiative and convective regions. Whether one in- 
cludes convective overshooting at the bound ary may also 
affect the depth of dredge-up. For example, IChiefB et al.l 
(|200l[ ) find that third dredge up does oc cur in their 
zero- m et allicity intermediate mass AGB stars. IChieffi et al.l 
treated t he convective boundaries according to the 
prescription of iHerwig et al.l ([l992), who use a mixing 
scheme so efficient that the composition discontinuity be- 
tween the two burning shells is smoothed out. This seems 
to indicate the efficiency of third dredge-up depends on 
the treatment of convection and the inclusion of extra- 



mixi ng mechanisms s uch a s convective overshooting. How- 
ever, ICil-Pons et al] l|200'H ) find that the total amount of 
mass dredged-up is very small even when overshooting is in- 
cluded. The STARS code employs an arithmetic mean when 
computing the mixing coefficient at the convective bound- 
aries. While this means we are mixing in material that lies 
outside the formal Schwarzschild boundary, it is a far smaller 
degree of overshooting than applied by other authors but a 
degree that leads to a stable composition and convection 
profile. We discuss the effects of overshooting in more detail 
in Section [l] We find that the inclusion of convective over- 
shooting in our model is enough to cause third dredge up to 

occur a t Z = 10~^ 

See lLau et all (|2008l ) for a more detailed discussion and 
description of the absence of third dredge-up. We notice that 
the 5Mq, Z = 10~^ star does not show third dredge-up. 
Also thermal pulses are growing weaker and weaker and 
eventually stop after 560 pulses in 5.5 x lO^j/r. Then the 
star enters a quiescent burning phase which lasts until car- 
bon ignites degenerately at the centre when the CO core 
mass is 1.36 Mq This agrees with the conclusion that su- 
pernovae Type 1.5 are likely to be the fate of these stars 
becau se of the faster c ore growth rate when thermal pulses 
cease (jLau et al.ll2008l l. 



3.4 2nd dredge up 

Fig. [2] shows the abundance of the CNO elements at the 
surface of each star at the beginning of the TP- AGB. We 
note that second dredge-up can lead to substantial enhance- 
ment of the CNO abundance. Of particular note is the 5 Mq , 
Z — 10~^ model, which has a surface CNO abundance that 
is about ten times higher than its initial metallicity. This en- 
hancement in metallicity is sufficient to allow the star to ex- 
perience a relatively normal TP- AGB. More CNO strength- 
ens the entropy barrier and decreases the H-tail from the 
H-burning shell. Note that at the same metallicity, the 3 Mq 
star experiences a CI phenomenon. Also, for the same mass, 
as we move to lower metallicity, the behaviour changes sud- 
denly as we drop the initial metallicity by just a factor of 
ten. 



3.5 Zero- metallicity models 

For stars at zero metallicity, the behaviour is similar to stars 
with Z = 10"®. In the 2Mq model, there is FDM at the first 
thermal pulse. For the 3 Mq model, the evolution is much 
more complicated. During the first pulse there is a CI phase 
but it is rather weak and does not increase the surface metal- 
licity by much. There is then FDM at the second thermal 
pulse. The 4M0 model encounters CI in the first pulse. The 
thermal pulses of the 5 Mq and 6 Mq model are not strong 
enough to cause any third dredge up and probably end their 
lives as type 1.5 supernovae (|Lau et al.ii2008i ). 



3.6 Comparison to models from literature 

ICampbell fc Lattanziol (|2008| ) have published a grid of mod- 
els for M/Mq = 0.85, 1.0, 2.0, 3.0 and the metallicity 
[Fc/H] = -6.5, -5.45, -4.0, -3.0 and Z = 0. Our re- 
sults for 2 Mq models agree well with their models. They 
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find FDM (referred as "Dual Sliell Flasli" by tfiem) at 
metallicity [Fe/H] = —4.0 or below but not at metallicity 
[Fe/H] = —3.0. This agrees with our result that the crit- 
ical metallicity for FDM to occur is between Z — 10~® 
to 10"^, corresponding to [Fe/H] between -4.3 and -3.3. For 
their 3M0 models FDM does not occur at [Fe/H] = -4.0 
or above but does when [Fe/H] = —5.45 and below. Our 
result is somewhat different. We only find the milder CI 
episode instead of FDM in all our 3 models except when 
Z = 0. CI begins to occur between Z = 10"^ and 10~^ 
([Fe/H] = —4.3 3.3) and this agrees with their metal- 
licity range for FDM at 3 Mq . For 1 Mq models we find 
the critical metallicity for fiash driven mixing to be between 

Z = IQ-'' and Z = 3 X 10"^ ([Fe/H] = -3.3 2.8). This 

agrees with the fact that FDM already occurs in their model 
with metallicity [Fe/H] = -3.0. 

One significant difference between the two sets of mod- 
els is the issue of hot bottom burning. Hot bottom burning 
does not occur for any of our models with initial mass 3 M© 
or below, while their 2 M© and 3 M0 models have strong 
hot bottom burning. Therefore, nitrogen is more abundant 
in their 2 M© and 3 M© models, while carbon is more abun- 
dant in our models. Hence, it is easier to form CEMP stars 
wit h [C/N] sa 1 from our models. 

iKomiva et al.l (|2007l ) give an up-to-date version of the 
dependence of FDM and the third dredge u p on the initial 
stellar mass and metallicity, formulated by iFuiimoto et al.l 
l|2000D . Our models show a few differences to their picture. 
Their critical metallicity for FDM to occur is [Fe/H] = —2.5. 
This is higher than the critical metallicity found in our grid 
of models ([Fe/H] = —3.3 to —2.8). Moreover, in their pic- 
ture, the critical metallicity for FDM to occur is the same 
for initial masses between 0.8 M© to 3.5 M©. This is not the 
case in our models (Fig. [TJ where the critical metallicity for 
FDM to occur in 2 M0 stars is lower than that of 1 M0 and 
1.5 Mq stars. Our mass range for FDM is also smaller than 
their picture. 

Another important difference is the nitrogen surface 
abundances. Their picture shows that extremely metal-poor 
stars between 3.5 Mq to 5M0 could be carbon rich but ni- 
trogen normal. However, most of our 4 — 6 M© models are 
nitrogen rich due to hot dredge-up. We agree that there are 
metal-poor stars that have no third dredge up but the mass 
range from our models is much smaller than theirs. 



4 THE EFFECT OF OVERSHOOTING 

The above behaviour raises an interesting question: what 
is the effect of the inclusion of convective overshooting on 
the behaviour of the models? To examine this, we have con- 
structed a grid of models with convective overshooting. We 
use the prescription of Schroder et al. ( 1997) and their over- 
shooting parameter, 5ov ~ 0.12. This overshooting is applied 
up to the beginning of the TP-AGB in all models made. 

Fig. [12] shows the abundance of the CNO elements at 
the surface of each star at the beginning of the TP-AGB. As 
with the non-overshooting case, we find that second dredge- 
up can lead to substantial enhancement of the CNO abun- 
dance. The effect is even more substantial with the inclusion 
of overshooting. This is what we would expect to find. We 
also note that we get a substantial elevation of the CNO 



abundance at lower initial masses. With overshooting our 
5M0, Z — 10~* model is significantly enhanced in CNO el- 
ements. The CNO abundances of 4-6M0, Z = 10"^ -10"* 
models with overshooting prescription is higher than models 
without overshooting by about one order of magnitude. 

Unfortunately we are not able to evolve some models 
along the whole TP-AGB with convective overshooting, ow- 
ing to as yet unresolved numerical difficulties during hot 
dredge-up for some stars. These models fail to converge at 
the onset of hot dredge-up if overshooting is applied. We 
therefore evolve without the overshooting just before the 
start of hot dredge-up. Fortunately we can evolve some mod- 
els, such as the 4M0, Z = 10~^, through a significant num- 
ber of thermal pulses with overshooting and, based on these 
models, we believe the inclusion of overshooting during the 
TP-AGB phase does not significantly affect the subsequent 
evolutionary behaviour. 

Models with overshooting have higher core masses when 
they enter the TP-AGB phase because of the more extended 
convective cores during core He-burning. Also, the amount of 
CNO elements dredged up is higher when overshooting is in- 
cluded. As a result, a model with a given mass and metallic- 
ity with overshooting behave more like a model with higher 
mass without overshooting. Comparing the evolutionary be- 
haviour with overshooting and without overshooting, we see 
that the major effect is to shift the boundary of different 
regimes of behaviour by mass but not metallicity (Fig. I13|) . 

We do not find any flash driven mixing or carbon in- 
gestion in stars with Z — 10"* as in the case of models 
without overshooting. At Z = 10~^ the overshooting mod- 
els of 2M0 show the occurrence of a relatively weak CI, 
which is absent in the models without overshooting. The 
overshooting prescription of STARS code mixes regions where 
Vr — Vad > —0.12. Therefore, regions that were close to in- 
stability can become convective. As a result, a small convec- 
tive region is opened up and carbon is mixed into proton-rich 
region. 

Unlike models without overshooting, models from Z — 
10"^ to Z = 10"'^ do not show any signature of FDM. It 
is possible that the maximum mass for a model with over- 
shooting to have FDM is just below 2 M0 for metallicities 
between 10~® and lO"'' because a strong CI episode is found 
in the 2M0 models at these two metallicities. Flash driven 
mixing does occur in our 2 Mq model with metallicity at 
10"^ 

There is a general trend that the upper mass limit for 
CI increases as metallicity decreases. At Z — 10~^ only the 
2M0 model shows a weak CI episode. CI occurs in both the 
2M0 and 3M0 models at Z = 10"® and in the 2M0, 3M0 
and 4M0 models at Z = 10"'^. At Z = 10"*, the 3M0 and 
4 M0 models still encounter CI but the 2 M0 model has a 
flash-driven mixing episode. 

It is worth noting that models of 6 — 7 M0 ignite carbon 
before any thermal pulses and become super- AGB stars. The 
5M0 models at Z = 10"** and Z = 10"^ also enter the 
SAGB phase before any thermal pulses and we do not model 
the evolution of these stars. This is in agreement with the 
general effect of overshooting which causes models to behave 
like higher mass models. 

We also flnd hot dredge-up in most of the 3 — 5 M0 
models, compared to the mass range of 4 — 6 M0 in the non- 
overshooting models. All of the 4M0 models encounter hot 



Low-Z AGB stars 11 




Mass/M0 


metallicity 


c 


14N 


16 o 


A 


3.0 


10-4 


2.4 X 10-4 


2.6 X 10-3 


3.0 X 10-4 


3.7 


3.0 


10-6 


2.6 X 10-4 


3.0 X 10-3 


2.8 X 10-4 


3.3 


4.0 


10-4 


7.1 X 10-4 


8.6 X 10-4 


8.1 X 10-8 


5.7 


4.0 


10-5 


3.7 X 10-8 


8.7 X 10-4 


7.4 X 10-8 


7.4 


4.0 


10-6 


3.8 X 10-8 


9.3 X 10-4 


9.8 X 10-8 


5.3 


4.0 
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5.0 
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5.0 
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1.0 X 10-8 


4.5 


5.0 
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2.7 X 10-8 


6.2 X 10-4 


5.8 X 10-8 
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Table 1. Surface abundances by mass fraction after the first hot dredge-up episode and dredge-up efficiency X of the hot dredge- up 
episode for models with overshooting. 



dredge up. While the two 5 M© models with Z = W~* and 
10^'' enter the SAGB branch, the other three 5 Mq models 
with lower metallicities also encounter third dredge up. For 
the five 3Mq models, we only find hot dredge up in the 
Z = 10"4. Z = 10"*^ models. Nevertheless, we believe hot 
dredge up should occur in all 3 M0 stars with metallicity 
ranging from 10^4 to 10^**. We did not find the occurrence 
in some of the SM© models because of the earlier break- 
down of the evolutionary code owing to unsolved numerical 
problems. 

Hot dredge up occurs at around the 5th pulse at metal- 
licity of 10~4 and 10~^. At low metallicity hot dredge up 
tends to occur in later pulses. In the SM©, Z = 10~* model 
it occurs at the 12th pulse. The dredge up efficiency A, 



defined as the mass dredged up divided by the core mass 
growth during intcrpulsc phase, is highest for 4M0 models 
in general (see table 1). The dredge up efficiency seems to re- 
main relatively constant between W~* and lO"'^ but is lower 
at Z = 10~*. For 5Mq models the dredge up efficiency is 
also the same for the Z = 10"^ and Z = 10"^ models but 
is lower in Z = 10~*. For the two SM© models, which show 
hot dredge up, A is also similar. In most cases, the code fails 
to converge at the next pulse after the hot dredge up but 
in a few, such as the 4Mq, Z = 10~^ model, we manage 
to evolve for two pulses with hot dredge up. The dredge up 
efficiency is much lower in the 2nd pulse, 1.5 compared to 
5.3 in the previous pulse. 

Also of interest is that the surface abundances after hot 
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Figure 13. The AGB behaviour of the models as a function of mass and metallicity with overshooting. The symbols are T — third 
drodge-up, H - hot third dredge-up, F - FDM, CI - CI, X - no TDUP and S - CI before the AGB (presumably leading to an S-AGB 
star). 



dredge do not depend much on metallicity. For 4M0 mod- 
els, the surface abundances of carbon are almost constant at 
3.8 X lO"'^ for Z = 10"^10"®,10"^ models (see Table 1). 
For nitrogen and oxygen, the surface abundances also hardly 
vary. M, Z — 10~* model there is more carbon enhancement 
than for the 4 M© models at lower metallicities while the ni- 
trogen and oxygen surface abundances are roughly the same 
for models with different metallicities. At Z — 10~* the sur- 
face abundances of carbon and nitrogen are about half those 
at higher metallicities while oxygen is down by nearly a fac- 
tor of ten. The CNO abundances for our SM©, Z = lO"** 
and Z = 10~^ stars are almost the same while the Z — 10~* 
model enriches its surface with more CNO elements by sev- 
eral times. The surface CNO abundances are also very close 
in the 3M0, ^ = 10"* and Z = lO"*^ models. This indi- 
cates that hot third dredge up wipes out any differences in 
the initial metallicity of the stars and that subsequently the 
surface abundances mainly depend on the initial mass. 

Amongst our overshooting models we do not find any 
that do not experience third dredge-up. One possible reason 



is that the second dredge up manages to increase the surface 
metallicity, in particular the CNO elements in the hydrogen 
shell, above the threshold at which the hydr ogen burning 
shell can be fully sustained by the CNO cycle l|Chieffi et al.l 
I2OOII ). 



5 CEMP STARS 

From our models, we can see that there are a few phenom- 
ena that can affect the composition of the material dredged 
up, namely flash-driven mixing, CI and hot dredge-up. If the 
material from the envelopes of these stars is accreted on to a 
binary companion it becomes carbon enhanced. In this sec- 
tion, we compare the surface abundances of our AGB star 
models results from our code to the observed abundances of 
different types of CEMP stars to check what types of CEMP 
stars can be formed by the binary transfer scenario. Ob- 
jects used in this section are picked from the Stellar Abun- 
dances of Galactic Archaeology (SAGA) Internet database 
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ISuda et alj |2008| ) and are chosen because their observed 
abundances fit well the abundances of our models. 

IStancliffe fc Glebbee^ (j200^ have discussed the forma- 
tion of CEMP stars by mass transfer from AGB stars in 
the mass range 1 Mq to 3.5 M0 with the effect of thermo- 
haline mixing in the secondary. One of the problems they 
note with the abundances from the AGB models is that 
they do not predict substantial enhancements of carbon and 
nitrogen at the same time except in a very narrow mass 
range. This problem is not solved by the evolution of higher 
mass AGB models in this work. Although carbon can be 
brought to the surface in third dredge up in these higher 
mass models, the presence of hot bottom burning destroys 
most of it. For example, in the 5Mq, Z = 10~* star with- 
out overshooting, [C/Fe] is roughly 1 and [N/Fe] is greater 
than 3 after hot third-dredge up and hot bottom burning. In 
the overshooting model, the surface abundances of carbon 
and nitrogen of the 4 Mq model are roughly equal. In this 
model, [C/Fe] is equal to 1 .6 and [N/Fe] is equal to 2.2 when 
[Fe/H] is about —2.4. This could possibly fit the observed 
abundances of CEMP stars with strong nitrogen enhance- 
ment, such as HE 1031-0020. The 3M0, Z = 10"" model 
gives a surface abundances of [C/Fe] — 1.1 and [N/Fe] — 2.7 
and may fit star such as HE 0400-2030. However, even with 
overshooting, the mass range in which we predict substantial 
enhancement of carbon and nitrogen is indeed narrow. 

At a metallicity oi Z — 10"*, we do not find any flash- 
driven mixing or carbon ingestion in any of our models. 
When we move down to the metallicity oi Z — 10~^ , an 
AGB star could have flash-driven mixing or carbon ingestion 
or hot dredge-up depending on its mass. More importantly, 
enhancement of carbon, nitrogen and other elements differs 
for these three mechanisms. For stars with FDM, the violent 
nature of the flash causes a very efficient mixing from the he- 
lium burning region to the hydrogen burning region. There- 
fore, when the convective envelope dredges into the region 
that is previously mixed by the FDM, significant amounts of 
the CNO elements are dredged to the surface. Carbon can 
be enhanced by more than 1,000 times (Fig. |3]). Nitrogen 
and oxygen are not enhanced as much and the C/N ratio 
approaches 100. For example, in the star of 1.5 Mq with 
Z = 10~^ subsequent thermal pulses further increase the 
carbon and oxygen surface abundance but not nitrogen be- 
cause the star is not hot enough for hot bottom burning. As 
a result, the [C/N] ratio can reach 1.5 or higher. Also, the 
[C/Fe] could be as high as 3. This could explain the forma- 
tion of HE 1430-0919 which has [C/Fe] = 2.7, [N/Fe] = 1.6 
and [Fe/H] = -3.01. 

For stars with a CI episode, we find that there is a range 
of carbon and nitrogen enhancement dependent on the mass 
of the stars. It is possible for stars with CI to achieve car- 
bon and nitrogen enhancement five times less than the stars 
with FDM with a similar [C/N] ratio. In some cases, such as 
the 3Mq, Z = 10~^ model with overshooting, the amount 
of carbon dredged up is lower and [C/Fe] only reaches 2.0 
and [N/Fe] only 0.9. This could well explain the formation 
of stars such as HE 2330-0555 and HE 0017-^0055. On the 
other hand, it is possible that CI could cause the nitrogen 
abundance to be similar to the carbon abundance as in the 
case of 3M0, Z = 10~^ without overshooting after the en- 
hancement by CI in the third pulse (see Fig. [5]). Whether CI 
preferentially enhances carbon or nitrogen depends on how 



efficiently the star converts carbon into nitrogen and hence 
depends on the meiss and metallicity of the stars. It has to 
be pointed out that, in the case of 3M0, Z = lO"'^ without 
overshooting, the subsequent third dredge up by later pulses 
eventually causes the surface carbon abundance to be higher 
than the surface nitrogen abundance. 

For higher-mass AGB stars at this metallicity, the hot 
dredge up could lead to simultaneous enhancement of carbon 
and nitrogen. However, while a significant amount of carbon 
is dredged up, most of it is converted to nitrogen by hot 
bottom burning, leading to a strong nitrogen enhancement. 
For example the 4Mq, Z = 10""" stars have surface abun- 
dances of [C/Fe] = 1.3, [N/Fe] = 3.1 and [Fe/H] = -3.3. 
Such nitrogen-rich stars have been searched for bu t there is 
an apparent dearth of them (| Johnson et al.ll2007l ). As the 
surface abundances of our models after hot dredge-up does 
not depend on metallicity, [N/Fe] increases with decreasing 
metallicity and stars with lower metallicity should have even 
stronger nitrogen enhancements if metal-poor AGB stars 
with hot dredge- up do exist. 

As we go down in metallicity to about 10~^, the second 
dredge up of high-mass AGB stars can reach deep enough to 
penetrate into regions where helium has been burned. Even 
though third dredge-up is weak or non-existent in such stars, 
carbon is already enhanced before entering the TP-AGB 
phase. For example, in the 5 Mq, Z = 10~* model with over- 
shooting, just after second dredge-up, carbon is enhanced 
by a factor of 10"^ while the C/N ratio is around 10. This 
could potentially be an important way of forming CEMP 
stars with extre mely low-metallicity, such as HE 0107-5240 
dLau et al.ll2007h . At the moment, we only observe two stars 
with an iron abundance of less than [Fe/H] < —5 (corre- 
sponding to Z — 10~^ if all the abundances were solar-scaled 
and the enhancement happened after these stars formed), 
HE 0107-5240 with [Fe/H] = -5.3 (|Christlicb et al. 20m 
and HE 1317-2326 with[Fe/H] = -5.45 (jFrebel et al.li2005h . 
However, this scenario cannot explain HE 1317-2326, be- 
cause of the similar observed enhancement of carbon and 
nitrogen. In fact, none of our models at such a low metallic- 
ity manage to produce a surface abundance with that C/N 
pattern. 

At such low metallicity, the carbon enhancement by 
high-mass AGB stars through second dredge up can be 
as significant as that of low-mass stars which go through 
FDM or CI. At metallicity of 10"^ and 10"^, the low-mass 
AGB stars have stronger carbon enhancement than high- 
mass AGB stars because FDM or CI leads to more efficient 
dredge up of carbon than third dredge-up alone. 

In conclusion, it is possible to explain CEMP stars with 
strong carbon enhancement and weaker nitrogen enhance- 
ment. They are most likely to be formed from low-mass AGB 
stars. At Z = 10~ these stars have third dredge up without 
hot bottom burning, so carbon is strongly enhanced. With 
metallicity around Z = 10~^ , both FDM and CI lead to sig- 
nificant carbon enhancement. It is harder to explain stars 
that have significant enhancement of both carbon and ni- 
trogen. While a narrow mass range of stars at Z = 10~* 
may be able to enrich carbon and nitrogen significantly 
through hot dredge-up, most models with hot dredge-up de- 
stroy too much carbon and enrich with too much nitrogen 
to fit the current observations. At a metallicity of 10~^ or 
below, nitrogen-enhanced stars are predicted by our models 
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but there is currently a dearth of observations of such stars. 
From our model, it is difficult to explain the existence of 
CEMP stars that have comparable [C/Fe] and [N/Fe] with 
[Fe/H] < —3. In order to produce stars with both carbon 
and nitrogen enriched from our model, we need hot bottom 
burning to be less efficient so less carbon is converted into 
nitrogen in high-mass AGB stars. If so, it may be possible 
to match those CEMP stars with both carbon and nitrogen 
enhanced. 

On the other hand, for some stars, there are brief 
times when both carbon and nitrogen are simultaneously 
enhanced, such as the 3M0, Z — 10~^ without overshoot- 
ing just after its CI episode but before hot dredge-up. Sim- 
ilarly, an AGB star can be first carbon-rich after the onset 
of third dredge up and then become nitrogen-rich when hot 
dredge up kicks in. If mass is transferred from the AGB pri- 
mary to the companion at suitable time, it is possible for the 
companion to end up both carbon and nitrogen enhanced. 
However, this involves a fine tuning of the period of the sys- 
tem and is unlikely to make significant amount of carbon 
and nitrogen enhanced metal-poor stars. 

Finally, we note that the composition of the transferred 
material may not necessarily be the same as that of an AGB 
primary which passes through its full TP phase. For exam- 
ple, if the star is in a close binary its life may be truncated if 
the star fills its Roche lobe or if the presence of a companion 
enhances its stellar wind. Therefore the binary properties of 
the system and the evolution of the surface abundances with 
time are parameters that could affect the observed compo- 
sition of the secondary. Moreover, the composition of the 
transferred material is also not necessarily the same as what 
will be observed on the secondary. Accreted material can 
become mixed with that of the secondary either through 
the deepening of the convective envelope when the star un- 
dergoes first dredge-up, or via thermohaline mixing imme- 
diately after accretion. The consequences of these processes 
for the s urface composition of the seco ndar y are discussed in 
detail bv lStancliffe fc Glebbeekl (|2008l ) and lStancliff3 (|2009l ) 



6 THE FEASIBILITY AND PROBLEMS OF 
DEDUCING THE INITIAL MASS 
FUNCTION 

"What is the initial mass function of metal-poor stars?" 
is an important question that has no conclusive answer 
at present. From the abundances of the observed carbon- 
enhanced metal-poor stars, one could theoretically deduce 
whether they were formed by binary mass transfer and if so, 
the initial mass of the primary stars. As a result, it ought to 
be possible to deduce the initial mass function of metal-poor 
stars based on the estimated mass of the primary stars. 

One important indicator of the primary mass of a par- 
ticular CEMP star is the nitrogen abundance, in particu- 
lar the [C/N] ratio, which ranges from —1 to around 1.3 in 
the observed CEMP stars. As mentioned above, a low-mass 
AGB star that enhances carbon through FDM or CI has a 
high [C/N] ratio unless hot dredge up is involved. The C/N 
ratio is usually to 2. On the other hand, for the high- 
mass AGB stars, because carbon is converted into nitrogen, 
the C/N ratio is much lower. The [C/N] ratio ranges from 
around -2 to about 0. 



Therefore, based on the observed [C/N] ratio, one might 
theoretically deduce whether the primary star is a high-mass 
or low-mass AGB stars. We need more data coming from 
ongoing surveys such as SDSS SEGUE, from which seven 
CEMP stars axe recently picked and discussed in lAoki et al.l 
(|2008l '). If enough CEMP stars are observed, one might de- 
duce the shape of the IMF of low-metallicity environment 
by noticing whether low-mass AGB are more numerous or 
vice versa. 

However we currently cannot match all the observed 
CEMP stars to initial primary masses and so a reliable IMF 
shape cannot be constructed. In particular, the high-mass 
AGB stars cannot fit some of the CEMP stars with strong 
enhancement in carbon and nitrogen and whether we should 
attribute all those type of stars to high-mass AGB stars 
would change the result. Though we can draw a general 
conclusion that a high fraction of CEMP stars with high 
C/N ratio implies more low-mass AGB primary and vice 
versa. 

On the other hand, our model predicts that both C/Fe 
and the fraction of CEMP stars increase at lower metal- 
licity. It would be interesting whether such a trend ap- 
pears when more metal-poor stars are observed, particularly 
CEMP stars with metallicity [Fe/H] < —4. If such a trend is 
not found observationally or the observed abundances, par- 
ticularly the C/N ratio, do not fit well with the surface abun- 
dances of AGB models, it could be that there is a significant 
shortcoming to current AGB modelling. Another possibility 
is that AGB stars did not form and the initial mass function 
peaked at the massive end below a critical metallicity. If this 
is the case, other formation channels are needed to explain 
those CEMP stars. 

6.1 The s-process elements 

We have not inclu ded full nucleosynthesis routines 
jStancliffe et al.ll2005l ). so abundances of light elements are 
not studied. This nucleosynthesis can also indicate whether 
s-processing has occurred. Without it we do not know if 
stars with FDM or CI mechanism have any s-process ele- 
ments. This could be important to deduce the mass of the 
primary because there are two groups of CEMP stars, one 
with s-process elements and one without s-process elements. 



7 CONCLUSIONS 

We have investigated the evolutionary properties of 
extremely-metal poor AGB stars with Z = 10"'* to Z — 
10~*. We show that for stars with metallicity oi Z = 10"'*, 
the AGB properties are similar to those of higher metallicity 
stars, except for the occurrence of hot third dredge up. How- 
ever, as metallicity decreases, flash-driven mixing or carbon 
ingestion occurs for low-mass AGB stars. The highest metal- 
licity models for which such phenomena occur is 10~^. We 
propose CI is a mild version of flash driven mixing. 

High-mass AGB stars with masses greater than about 
4 Mq can undergo hot dredge-up. As a result, the surface 
C/N ratio of these stars is lower. Surface abundances after 
hot dredge up depend on mass but not metallicity. The criti- 
cal metallicity below which third dredge up does not occur is 
Z — 10~*. It is possible for some stars at this metallicity to 
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completely cease thermally pulsing later i n their evolutio n 
and end their fate as supernovae type 1.5. (|Lau et al.ll2008h . 
One model, 5M0, Z = 10~* without overshooting, is found 
to have this type of behaviour. However it does not occur 
with overshooting for stars with Z — 10~*. Overshooting 
also causes the models to behave like a higher mass models 
without overshooting. In general, the mass boundaries of the 
different behaviours are shifted by about 1 Mq . 

We show that the C/N ratio of low- mass AGB stars 
can be as high as several hundred while the C/N ratio of 
high-mass AGB stars can be less than one. The variation 
of carbon and nitrogen enhancement of CEMP stars could 
then be due to the different initial mass of the primary. 
However our models here can only explain the CEMP stars 
with C/N much larger unity and those with C/N around 
1/10. They struggle to explain CEMP stars with similar 
enhancements of carbon and nitrogen. Hot third-dredge up 
is needed to convert carbon to nitrogen but it is too efficient 
in our models, destroys too much carbon and produces too 
much nitrogen to fit the observations. 

Theoretically, based on the observed frequency of differ- 
ent C/N ratio of CEMP stars at low metallicity, it is possible 
to deduce the shape of initial mass function of AGB stars, 
as stars with large C /N ratio are more likely to accrete mass 
from low-mass AGB stars. 
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